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^! Abstract 

< 

_i. ' The strong energy dependence of the s-wave r]N scattering amphtude at 

CN ■ and below threshold, as evident in coupled-channels i^-matrix fits and chiral 

models that incorporate the Su A^*(1535) resonance, is included self consis- 
tently in 77-nuclear bound state calculations. This approach, applied recently 
in calculations of kaonic atoms and ^-nuclear bound states, is found to im- 
pose stronger constraints than ever on the onset of r^-nuclear binding, with a 
pj ■ minimum value of Re Qr^N ~ 0.9 fm required to accommodate an r^-'^He bound 

state. Binding energies and widths of ?7-nuclear states are calculated within 
several underlying r]N models for nuclei across the periodic table, including 
^^Mg for which some evidence was proposed in a recent COSY experiment. 
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1. Introduction 

/\ . Searches for meson-nuclear bound states have focused on K~ and r] 

c^ \ mesons, motivated by a general theoretical consensus that the near-threshold 

^A^ and r]N attraction generated by the s-wave resonances A(1405) and 
A^*(1535), respectively, translates into sufficiently attractive A^'-nucleus and 
?7-nucleus interactions. A corollary of this resonance dominance is a strong 
energy dependence of the underlying ^A^ and r]N interactions. Here we 
apply the lessons gained by handling the strong energy dependence of the 
near-threshold KN interaction in A'~-nuclear calculations [l| to ry-nuclear 
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bound-state calculations. Early calculations by Haider and Liu |2|, y] pre- 
dicted r^-nuclear bound states beginning with nuclear mass number A ~ 12. 
In these, as well as in a follow-up calculation J4|, a fairly weak rjN attrac- 
tion input was used, with Re a^Tv ^ 0.3 fm, where ar^N is the rjN scattering 
length. Several versions of coupled-channels chiral models give similar val- 
ues as well [a, la, LZI] whereas other models, particularly those using i^-matrix 
methods to fit ttN and 7A^ reaction data in the A^*(1535) resonance re- 
gion, e.g. [sl, y, |lO|, yield considerably stronger rjN attraction with values 
of Re ttriN ~ 1 fmo This might suggest that the onset of ?7-nuclear binding 
occurs already in the He isotopes for which strong final-state interaction pre- 



cursors have been noted in deuteron-initiated rj production [ll|, ll2|. A robust 



pattern of ry-nuclear bound states could yield useful information on the size 



of SU(3) flavor 77 — 77' mixing and about axial U(l) dynamics 13 
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Figure 1: Energy dependence of the s-wave rjN cm scattering amplitude friN{\/s), with 
scattering length o^at = frjN{^/sth) = 0.91(6) + i0.27(2) fm, as calculated by Green and 
Wycech (GW) fitting pion reaction and photoproduction data over a wide range of cm 
energies Ecm (=v^)- The spike of Re / is located at the t]N threshold, whereas Im / 
peaks close to the A^*(1535) mass. Figure adapted from Ref. [8|. 



Regardless of the strong model dependence of Re a^^r, all studies of the 
TjN system near threshold, y^Sth = ^^n + ^^77 ~ 1487 MeV, agree that both 
real and imaginary parts of the s-wave center-of-mass (cm) scattering am- 



^As for Im arjN, it is constrained by ttN — > rjN cross-section measurements, with values 
of Im a^jv ^ 0.2-0.3 fm in most theoretical analyses. 



Table 1: rjN scattering length a^^v (in fm) in three coupled-channels models used in the 
present work. Ml and M2 correspond to versions I and II, respectively, of rjN amplitudes 
from the recent chiral-model work by Mai et al. [7[, and GW denotes the if-matrix rjN 
amplitude due to Green and Wycech ^ shown in Fig. [T] 



Model 


Re CL-ql^l 


Im ariN 


Ml 


0.22 


0.24 


M2 


0.38 


0.20 


GW 


0.96 


0.26 



plitude fj^N decrease steeply in going subthreshold, as illustrated in Fig. [H 
Since the in-medium rjN interaction relevant to the evaluation of ?7-nuclear 
bound states involves subthreshold rjN configurations, a procedure for go- 
ing subthreshold is mandatory. Previous calculations focused on shifting the 
energy variable of /,,Af(^/i) or its in-medium version by a fixed amount be- 
low threshold: 5y/s = —30 MeV was found in Ref. |j| to provide a good 
approximation to a variety of off-shell effects, whereas 6^/s = —B^i, with 



B^i the r^-nuclear binding energy, was used in Ref. IJ]. The latter proce- 
dure requires a self-consistent calculation to ensure that the S^ argument of 
the input /,,Ar coincides with the -B^ output of the binding energy calcula- 
tion. However, it was shown in our recent studies of ii"~-nuclear dynamics 



15L lid . llTl . llSL Il9l . |20| that a more involved self-consistent calculation is re- 
quired to correctly implement the subthreshold energy dependence, and it 
is this self-consistent procedure that is applied here to calculate ?7-nuclear 
bound states. This procedure results in imposing stronger constraints than 
ever on the onset of 77-nuclear binding. 

Below we proceed to describe briefly the self-consistent procedure used to 
handle the subthreshold energy dependence of the rjN amplitude for bound 
nucleons, and its embeding into a dynamical Relativistic Mean Field (RMF) 
scheme which allows for the flrst time to consider the polarization of the 
core nucleus by the bound t] meson. To span a broad range of bound-state 
scenarios we apply our methodology to three distinct rjN amplitude models, 
with threshold values listed in Table [1] These amplitudes differ primarily 
in the value of the real part, while their shape below threshold exhibits 
a substantial decrease particularly for Im a^Ar, as illustrated for the GW 
amplitude model |8[ in Fig. [H We have calculated 77-nuclear bound states 



across the periodic table for these three amphtude models, as reported and 
discussed here for Is^ states. Finally, we also confront our results with a 



recent experimental suggestion of a ^Mg bound state [21 . 



2. Methodology 



In close analogy to the latest calculation of K~ nuclear bound states [18 
we calculate t] nuclear bound states by solving self consistently the Klein- 
Gordon (KG) equation 

[ v^ + ul -ml- n,K, p) ] v^ = 0, (1) 

where u^^ = Un — ir^/2 and w^ = m^ — B^, with B^ and F^ the binding 
energy and the width of the r] nuclear bound state. The self-energy operator 
n^ is related to a density- and energy-dependent optical potential K7 which 
is given by the following "tp" form: 

11^(0;^, p) = 2ujnVr, = -AnFriNiVs, p)p, (2) 

where s = {E^ + E]^Y — (j)^ + p^Y is the Lorentz invariant Mandelstam 
variable s which reduces to the square of the total r]N energy in the two- 
body cm frame and F^at is the in-medium rjN s-wave scattering amplitude 
in the lab system. Note that for A^ 1 the lab system approximates well the 
?7-nucleus cm system. Our in-medium F^jv accounts for Pauli correlations in 



the Ericson-Ericson multiple-scattering approach, as reformulated in Ref. [22 



and used recently in Ref. [20 



FM^s,p) = ^-^^^^ , e(p) = ^, (3) 

i + e(p)V(v^)p ^pV ^ 

where fr,N{,y/s) = {y/s/mN)frjN{y/s), with the kinematical factor y/s/mj^ 
transforming / from the two-body cm frame to the lab /, and where pp is 
the local Fermi momentum corresponding to density p = 2p|./(37r^). Note 
that Fj^N{y/s,p) — i- fnN{y/s) upon p — )• 0, as required by the low-density 
limit. Extensions of Eq. ((3)) to coupled channels and inclusion of self energies 
are discussed elsewhere [23.] . 

The in-medium y/s = -^/sth + S^/s energy argument of F^at in Eq. ([3]) is 



density dependent, with a form used recently in K -atom studies [20 



6,/^ ^-Bn-- ^nB^^ - ^NTN{-f'^ + ^r^Re K,(5^, p), (4) 

P Po Po 



where ^Ar(r;) = f^N {■q) I ij^N + f^ri) , B]\r ^ 8.5 MeV is the average nucleon bind- 
ing energy, T^r ^ 23 MeV is the average bound-nucleon kinetic energy in the 
Fermi Gas Model, and po (p) is the maximal (average) nuclear density. The 
Vrj term in expression (jl]) originates from the nonzero in-medium momentum 
(squared) of the r] meson in finite nuclei, and is related in the local density 
approximation to the r^-meson kinetic energy [15]. This contrasts with the 
common assumption Pr^ = made in nuclear matter calculations. The depen- 
dence of Vrj on energy through 5^ and on density p is explicitly marked in this 
expression. Eq. (jl]) satisfies the low-density limit requirement that 6^/s — )■ 
with p — )■ 0. Note that for an attractive K; and as long as p 7^ 0, the shift 
of the two-body energy away from threshold is negative definite, 6y/s < 0, 
even as i?^ — )■ 0. For a given 5^, neither a/s nor Vr^ can be evaluated sepa- 
rately, implying that V^ is to be constructed self consistently together with 
y/s, which takes typically about 5 cycles of iteration. Once V^(i?j,,p) has 
been determined, it is used in the KG Eq. ([1]) to solve for the binding en- 
ergy eigenvalue Bi^ in the 77-nuclear n(. single-particle state. While varying 
the value of B^ in this process, the self-consistent requirement Eq. ^ is 
imposed at each step of the calculation of the eigenvalue. 
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Figure 2: Subthreshold rjN energies probed by the rj nuclear potential as a function of 
the relative nuclear RMF density in Ca. Each of the three curves was calculated self- 
consistently according to Eq. ^ for a specific version of rjN subthreshold amplitude 
model, see text. 



In Fig. [2] we show the downward subthreshold energy shift S^/s = E — Et^ 
as a function of the nuclear RMF density p in Ca, calculated self consistently 
according to Eq. (jl]) for riN amplitude models Ml, M2 and GW (see Ta- 
ble [1]). The hierarchy of the three curves reflects the strength of the input 
Re fr)N{.\fs) in the subthreshold region, with threshold values listed in Ta- 
ble dJ It is clear that downward energy shifts of up to ~55 MeV are involved 
in the present self-consistent calculations. 

3. Results and discussion 

The methodology described in the last section was used to solve the KG 
equation ([1]) for r^-nuclear bound states across the periodic table. In this 
Letter we highlight the systematics of the Is^ bound state and compare our 
treatment of subthreshold energy dependence with previous studies. A more 



detailed account plus extensions are given elsewhere 23j. Three represen- 
tative rjN amplitude models Ml, M2 and GW (see Table [1]) are employed 
here in order to span a wide range of tjN interaction strengths. Our main 
results are shown in Fig. [3] for binding energies i?^ and widths F^^ calculated 
for Is^ nuclear states in core nuclei from ^^C to ^°'^Pb. RMF equations of 



motion, along with the KG equation ([T]), are solved self consistently [24 



thereby allowing for core polarization by the rj meson (see Ref. [16| for the 
latest application to K~ mesons). The core polarization effect on 5^ and F^ 
was found in all cases displayed here to be less than 1 MeV. Therefore, the 
use of static nuclear densities is acceptable for not-too-light nuclear cores. 

Inspection of the l.h.s. of Fig. |3] reveals that for each of the three input 
rjN amplitude models the binding energy increases with A and appears to 
saturate for large values of A. As in Fig. [2l here too the hierarchy of the 
three curves reflects the strength of the Re frjNiVs) input in the subthreshold 
region, with threshold values listed in Table [H The Ml and M2 amplitudes 
are too weak to produce a Is,, bound state in ^^C, with the onset of binding 
for the weaker Ml amplitude deferred to around ^°Ca. Of our three represen- 
tative amplitudes. Ml is the closest one on shell to the Haider-Liu standard 
amplitude [2|, y, \^ which was used by these authors to argue for ^^C as the 
approximate onset of r] nuclear binding. In contrast. Re fr,N{,y/s) of the GW 
model is sufficiently strong to bind the Is^ state in ^^C and in lighter core 
nuclei, in spite of the suppression it undergoes here by forming its in-medium 
version and dealing with its energy dependence. The GW amplitude model 
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Figure 3: Binding energies (left) and widths (right) of Is,, nuclear states across the periodic 
table calculated self consistently using the Ml, M2 and GW subthreshold ijN scattering 
amplitudes within a dynamical RMF scheme, see text. 



even admits a Isrj bound state in ^He with as low a binding energy as 1.2 MeV 
and a width of 2.3 MeV, both calculated using a static ^He density. 

Inspection of the r.h.s. of Fig. [3] reveals a trend for the three curves of 
calculated widths which is opposite to that observed on the l.h.s. for the 
calculated binding energies. Here, the GW model produces relatively small 
widths of order 4 MeV uniformly across the periodic table, whereas Ml and 
M2 give larger widths, particularly Ml with widths of order 20 MeV. This 
reflects partly the energy dependence of Im friN{y/s) in the subthreshold 
region, which is quite distinct in each one of the three amplitude models, 
and partly the difference in the in-medium renormalization arising from the 
Re friN{\/s) input. The latter point is readily understood by noting in Fig. [2] 
that the largest values of subthreshold downward energy shift are due to the 
GW subthreshold amplitude. This causes a particularly large reduction in 
the strength of the Im friN{y/s) input for the GW amplitude model. 

Focusing on a given core nucleus, we show in Table |2] results of static- 
density calculations in models Ml, M2 and GW of the Is,, state in ^^Mg with 
(Yes) and without (No) employing the in-medium modification of Eq. (|3]). 
The self-consistency requirement imposed by Eq. (jlj) is used and comparison 
is made with two other procedures applied in previous studies to incorporate 
energy dependence. These are (i) a fixed 30 MeV downward shift applied to 
the free-space r]N amplitude by Haider and Liu (HL) [4]; and (ii) shifting 



Table 2: Static calculations of Is,, binding energy (B) and width (F) in ^^Mg, using three 
rjN amplitude models (Ml, M2, GW) with (Yes) and without (No) medium corrections 
from Eq. ([3]), for several procedures of treating the energy dependence of fr^N- Energies 
and widths are given in MeV. 

Eq. §i> subthreshold Ml M2 GW 







B, 


r,, B^ 


r. 


B, 


r. 


No 


5^ = -30 


- 


- 3.0 


11.2 


31.2 


10.0 


No 


(5v^ Eq. (iD 


- 


- 3.2 


10.6 


23.8 


7.4 


Yes 


6^/s = -B^ 


3.8 


21.7 8.3 


13.2 


19.4 


5.8 


Yes 


5^s Eq. (H 


- 


- 2.5 


7.4 


14.8 


3.9 



down self-consistently the energy argument of the in-medium f^^^ arnplitude 
by the resultant i?^, as implemented by Garcia- Redo et al. (GR) [ij]. In 
the first two lines of the table we compare the HL procedure with ours, 
using free-space amplitudes. Both do not produce binding for the weakest 
amplitude Ml and practically agree for M2, while disagreeing significantly 
for the strongest GW amplitude. In the last two lines we compare the GR 
procedure with ours, using in-medium amplitudes. The GR procedure is 
found to give significantly higher binding energies and widths than ours for 
all amplitude models tested here, particularly for the weaker M amplitudes 
where it is the only one that produces a Is,, bound state for Ml. 

Of the three models used by us with in-medium amplitudes in Table [2] 
(last line), only GW provides B^, which is comparable with 

S^"P(2jMg) = 13.1 ± 1.6 MeV, r"P(2jMg) = 10.2 ± 3.0 MeV, (5) 

deduced from the following ^^Mg interpretation of a peak observed by the 
COSY-GEM Collaboration ^■. 

p + 27 Al ^ 25^g ^ 3gg 

^ (tt- + p) + X, (6) 

with a decay induced hj r] + n ^ n^ + p. Hence, if this peak assignment to 
a Isri state is correctp then the underlying threshold value Re ajjN must be 



^This has been contested recently by Haider and Liu who offered a different interpre- 
tation of the observed peak [25| . 

8 



rather large, close to 1 fm. Other procedures listed in Table [2] for treating 
the subthreshold r]N energy dependence require considerably smaller values 
of Re a^Tv- Finally, the relatively small value of width F produced in the 
GW model should not be viewed as too restrictive since the total width must 
be larger than given in these models, owing to true rjNN absorption and 
two-pion production rjN — )■ -r-kN processes that are not accounted for by the 
models considered in the present work. 

4. Conclusions 

In this work we have demonstrated the importance of, as well as the sub- 
tleties involved in constructing self-consistent 77-nucleus optical potentials 
that incorporate the strong subthreshold energy dependence of the under- 
lying rjN scattering amplitude. Of the three rjN amplitude models studied 
here self-consistently, even the relatively weak attraction in model Ml with a 
threshold value Re a^A? ~ 0.2 fm requires going down to about 35 MeV below 
threshold, as shown in Fig. [21 in order to calculate reliably the ry-nuclear op- 
tical potential Ky(p) at central nuclear densities. This downward energy shift 
exceeds by far the downward shifts —B^, with i?^ < 20 MeV encountered 



in the self- consistent calculations of Garcia- Redo et al. [ij]. The relatively 
large downward energy shifts in the present approach together with the rapid 
decrease of the free-space and in-medium rjN amplitudes lead to smaller than 
ever binding energies calculated in comparable models ^ [ij]. Thus, ^^C 
bound states are unlikely in models with threshold values Re a^^v ^ 0.5 fm, 
and as large a value as Re ar^N ~ 0.9 fm is required to reproduce the ^^Mg 



bound-state hint from the recent COSY-GEM experiment [21 

A value of Re a^^v ~ 0.9 fm is likely to yield a near-threshold ^He bound 
state, as found by us using the GW amplitude model, but it is short of binding 
^He. Stretching the limits of optical potential usage down to these light 
systems is of course questionable, and corresponding few-body calculations 
are highly needed to resolve such issues. Nevertheless, if one applies our 
subthreshold self-consistency scheme to r^-^He low-energy scattering, then 
a very large imaginary part that might indicate a nearby virtual state is 
found for the ?7-^He cm scattering amplitude when using the GW amplitude 
model. This large imaginary part might be associated with the strong final- 
state interaction effects observed for the ry-^He system ll|. In contrast to 



previous estimates that assigned a value of Re ar/N ~ 0.5 fm to describe such 



occurrence [26|, in our self-consistent calculations it requires substantially 
larger values, more likely around 0.9 fm. 
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